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Abstract 

In this work we demonstrate that the Yang-Baxter algebra can also be employed in 
order to derive a functional relation for the partition function of the six vertex model 
with domain wall boundary conditions. The homogeneous limit is studied for small 
lattices and the properties determining the partition function are also discussed. 
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1 Introduction 

The six vertex model with domain wall boundary conditions was introduced in [1] and since 
then many connections with enumerative combinatorics and orthogonal polynomials theory 
have been unveiled [2]. For instance, the problem of alternating sign matrices (ASM) and 
domino tilings enumeration is known to have a close relationship with the six vertex model 
with this particular boundary condition (3j HJ El E] . Moreover, the partition function of this 
model was also shown to correspond to a Schubert polynomial in [7] and to a KP r function 
in g]. 

From the physical point of view the study of the partition function of this model reveals 
us an interesting phenomena. Using the determinant representation found by Izergin [H] and 
a Toda chain differential equation [TU] satisfied by this partition function, it was shown in 
[TTj that the bulk free energy of the six vertex model with domain wall boundaries in the 
thermodynamical limit differs from the one with periodic boundary conditions. This fact 
has been also discussed in [12J and raises the issue of the sensitivity of the six vertex model 
in the thermodynamical limit with boundary conditions. 



Domain wall boundary conditions were first introduced and studied within the scope 
of the Quantum Inverse Scattering Method P] , and they emerge naturally in the calculation 
of correlations functions of quantum integrable systems [13]. On the other hand, integrable 
systems with an underlying Yang-Baxter symmetry have also been tackled by functional 
equations methods, which are intimately connected with Baxter's commuting transfer matrix 
approach [TJ]. However, since functional methods usually do not provide the eigenvectors of 
the system, it is usually claimed that the calculation of correlation functions is out of reach 
for functional equations methods. 

The use of the Yang-Baxter algebra in order to obtain functional equations was first 
introduced in [15] and the aim of this paper is to show that the Yang-Baxter equation can 
still be explored in order to obtain a functional relation for the partition function of the six 
vertex model with domain wall boundaries. 

This paper is organized as follows. In the section 2 we describe the six vertex model with 
domain wall boundaries and its construction in terms of the Yang-Baxter algebra elements. 
In the section 3 we derive a functional equation determining the partition function of the 
model and in section 4 and 5 we study its homogeneous limit and some particular solutions. 
Concluding remarks and open questions are discussed in the section 6 and in the appendix 
A and B we present some extra results and technical details. 



2 The six vertex model with domain wall boundaries 

Vertex models in Statistical Mechanics were first introduced by L. Pauling aiming to describe 
the residual entropy of ice [16] . These models are described in terms of a matrix £ containing 
the statistical weights of the possible vertex configurations. 

In a two dimensional rectangular lattice consisting of M horizontal and L vertical lines 
we have Lx M intersection points and the intersection point of the i-th horizontal and the j- 
th vertical line together with the four connecting edges is referred as a vertex, as represented 
in the Fig. 1. 

The edge variables {ajj, a^+i, fyj, fy+ij} characterizes the configuration of the vertex 
at position while C^' 3+ ^. denotes its statistical weight. The partition function of the 
model is thus obtained by summing the product of the statistical weights over all possible 
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Figure 1: (a) Rectangular lattice with M lines and L columns, (b) The vertex at position with its 

corresponding statistical weight. 



configurations, i.e. 

M L 

z znnc ; w 

{<*,/?} i=l j=l 

For a more detailed discussion about vertex models and their applications in Statistical 
Mechanics see for instance [17] and references therein. 

In order to evaluate the sum ([T]) it is usually assumed a particular boundary condi- 
tion and the concept of domain wall boundary conditions for the inhomogeneous six vertex 
model was introduced by Korepin in the Ref. pQ. Turning our attention to the six vertex 
model, each edge variable and can assume two possible configurations under a certain 
restriction known as ice rule. These assumptions result in the following £-matrix 



£(A) 



a(A) 






\ 



(2) 





b{\) c(A) 
c(A) 6(A) 
V U a(A) I 

containing the statistical weights a, b and c of the allowed configurations. The weights are 
explicitly given by a(A) = sinh(A + 7), b(\) = sinh (A) and c(A) = sinh (7) where 7 is the 
anisotropy parameter and the complex variable A parametrizes the Yang-Baxter integrable 
manifold 

n 2 _|_ h 2 _ r 2 

A A = cosh (7). (3) 
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Strictly speaking, the £-matrix (j2]) satisfies the Yang-Baxter relation 

£ 12 (A - A*)£i3(A - iy)C 23 ^ - v) = £ 2 3(> - z^)£i 3 (A - i/)£i 2 (A - jx) (4) 
where £jj G End (V^ £g> Vj) and Vi = C 2 . Consequently, the monodromy matrix 

T(A, {/i fc }) = £^i(A - ^i)C A2 (X - fa) ■ ■ • £al(A - n£) (5) 
satisfies the relation 

r(x - u) r(A, {fi k }) ® 7>, {M) = 7>, K» ® t(a, {^» #(a - ^) (6) 

where R(X) = P£(A) and P is the standard permutation matrix. The relation (Q is com- 
monly referred as Yang-Baxter algebra and, together with the relation (j3J), they constitute 
the basis of the Quantum Inverse Scattering Method [T8J Q2] . 

Boundary conditions are an important ingredient in the formulation of lattice models 
and the case with domain wall boundary conditions introduced in pQ requires that the 
boundary edges have a particular configuration respecting the six vertex model symmetry. 
For the six vertex model each edge variable {otij, fiij} can assume two possible configurations 
which can be conveniently denoted by arrows pointing inwards or outwards, i.e. 

ctij = — > or ^~ 
= I or | 

In this way domain wall boundary conditions consist of the restrictions 

Oii,l = -> &i,L+l = <- 

0i,j = I Pl+ij = t (8) 

where now we are considering a square lattice with M = L. In order to make this situation 
more clear and intuitive, we have depicted a possible six vertex lattice configuration with 
domain wall boundaries in the Fig. 2. 

For the six vertex model the monodromy matrix (j5J consist of a 2 x 2 matrix with 
operator valued entries. Within the framework of the Quantum Inverse Scattering Method 
it can be conveniently denoted as 

r ( A,w)=h A ' K}) B(XM) ) ( 9 ) 

\C(X,{H k }) D(X,{fi k })J 
4 



(7) 




Figure 2: A configuration with domain wall boundary conditions. 



and it was shown in [T] that the partition function (pQ) of the inhomogeneous six vertex model 
with domain wall boundary conditions can be written as 

L 



Z DW = (Q\l[B{X j ,{fjL k }) |0) 

i=i 

or equivalently 

L 

Z DW = (0\l[C(X„{fi k })\0) . 
where |0) and |0) denote the usual ferromagnetic states 

io>=|(;) ^ ,o)^|(; 



(10) 



(11) 



(12) 



In the next sections we shall demonstrate how the Yang-Baxter algebra can be employed in 
order to obtain a functional relation for the partition function ffTUj) . 



3 Yang-Baxter algebra and functional relations 

The Yang-Baxter algebra is a corner stone of the Quantum Inverse Scattering Method and it 
has been successfully explored in order to construct integrable systems and to extract exact 
results from them. More recently, it was shown in [15] that the Yang-Baxter algebra can 
also render functional relations determining the spectrum of spin chains with non-diagonal 
twisted and open boundary conditions. 

In the Ref. pQ Korepin obtained a recursion relation for the partition function (fTUl) 
whose solution was later on given by Izergin [19] in terms of a determinant formula. Our 
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purpose here is to demonstrate that the Yang-Baxter algebra also provide us a functional 
equation determining the partition function of the six vertex model with domain wall bound- 
aries. In order to do so we first need to recall some properties exhibited by the six vertex 
model monodromy matrix 

M M 

Let |0) • and |0) ■ be the states and I respectively, acting on the j-th space of 

W W 

the tensor product V\ <8> V2 <E> • • • <8> Vj <g> • • • <8> Vl- The action of £4, on these states yields us 
a triangular matrix. For instance we have 

/ a(A) 10), f \ 
Las A \% =\ K)U ' (13) 
V & ( A )I W 

and analogously 

/ 6(A) 10), \ 
4* A , = ^ ' h , (14) 

V t «(A)|0) iy / 

where the symbol f stands for a non-null value. 

Due to its definition, the monodromy matrix (jSJ) inherits the triangular properties f lT5|) 

and ( FT4|) . In this way the monodromy matrix elements satisfy the relations 

L 



A(X, {fi k }) |0> = J] a(X - n s ) |0> B(X, {„»}) |0) =f 



i=i 



L 



C(X, {&}) |0> = D(X, fa}) |0) = J] b(X - ^) |0> , (15) 

i=i 

and 

L 

A(A, {/**}) |0) = II fe ( A - ^) |0> B(A, W) |0> =0 

L 

C(A, {/**}) |0> = f 0(A, {M) |0> = II a ( A - Mi) 1°) • ( 16 ) 

i=i 

Another fundamental ingredient in the method considered here is the Yang-Baxter alge- 
bra. The relation fl6]) encodes several commutation rules for the elements of the monodromy 
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matrix and among them we shall make use of the following ones 
A(X,{fi k })B(u,{fji k }) : 

D(A,{M)#MM) 

[C(\,{fM k }),B{v,{n k })] 
B(\,{fi k })B(u,{vL k }) 



B(u, { Mfe })A(A, fa}) - C ^—^B(X, {^ k })A{v, fa}) 

B(u, {^ k })D(X, {/i fe }) - ^—^±B(\, {^})D(u, {/i k }) 

[A(u, {fi k })D(X, {/x fe }) - A(X, M)D(u, {fi k })} 
B(u,{^ k })B(X,{fi k }) . (17) 



a[y — 


A) 


h{y- 


A) 


a(X — 


v) 


b(X- 


v) 


c(A- 


v) 


b(X- 


v) 



In the Refs. [18J and [lj it was shown that the operator B(X, {^k}) plays the role of 
raising operator with respect to the pseudo- vacuum state |0) while the operator C(A, {fJ> k }) 
acts as a lowering operator. Furthermore, it was also demonstrated in [TJ that 

n n n 

C(X ,{fi k })l[B(X t , {//*}) 10) = ^M»[]%,W)|0) 



with 



i=l 



i=l 3=1 



+ yi ^(A ,{^})n 5 ( A ^{^})i°) 

l<i<j<n 



1=1 



b(Xi - X 
c(A - Aj) 



i=i 

L 



11 b(Xi - A fe ) b(X k - A ) 

k— 1 

a(A - A fc ) a(X k - \) 



k^i 
n 



k=x 
k+i 



b(X - A fe ) b(X k - Xi) 



(19) 



N 



c(A - Xj) c(Xi - A ) a(Xj - Xi) 



a(Xj — X m ) a(X m — Xi) 



b(X - A,) b(X t - A ) b(Xj - A ? .) J! a ( A * ^ II ^ _ Am ) b(X m - A,) 



c(A - Xj) c(Xj - A ) a[Xj - Xj) 
b(X - A,) b(Xj - A ) b{\ - Xj) 



Yl a (Aj - Vi)b{Xi - hi) ] [ 



a (A, 



A m ) a(A m — Xj) 



i=i 



\ b(Xi - X m ) b(X m - Xj) 

m^i,j 

(20) 



for any number n of operators B(Xi, {/J. k })- The demonstration of (|18|) - (|20p only makes use 
of the commutation rules (IT7|) together with the relations f[To"j) . 

At this stage we have gathered most of the ingredients required to obtain a functional 
relation for the partition function (flOj) . In order to proceed we look to the relation (Fl8|) with 
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n = L + 1 and act with the dual vector (0| on its left hand side. By doing so we obtain the 
following relation 

L+1 L+1 

(0| C(A , W) J] B(Xi, {fi k }) |0) = M i z ( x ^ ■ ■ ■ > A -i' Ai+i, . . . , A L+1 ) 

1=1 i=l 

+ -^(Ao, Ai, . . . , Aj_i, Aj + i, . . . , Xj-i, Aj + i, . . . , Xl+i) 

i<*<i<L+i 

(21) 

where Z(Xi, . . . , A/,) denotes the partition function (TTUj) omitting the dependence with the 
variables {/ifc}. 

n 

We shall refer to Y\ B(Xi, {fJ-k}) |0) as Bethe vectors and it was also shown in [TJ that 
i=i 

L 

[I%fc})|0}-2(A lr .,A L )|0> (22) 
i=i 

which implies in 

L+1 

IJfl(A i ,{A»*})|0>=0 (23) 

3=1 

due to the relations f JT6|) . The Bethe vectors enjoy the property of being highest weight 



SU(2) vectors [2D1 EE] and a careful examination of f[2"Tj) . taking into account the relation 
(l23p . reveals that its left hand side vanishes. In other words, the consistency of the Yang- 
Baxter algebra with the highest weight property of the Bethe vectors results in the following 
functional equation for (flOj) . 

L+1 

Mj Z(X\, . . . , Aj_i, Aj + i, . . . , A^ + i) 

i=i 

+ Nji Z(Xo, Ai, . . . , Aj_i, Aj + i, . . . , Xj-i, Xj+i, . . . , Al + i) = , (24) 

l<«<j<L+l 

where the functions Mi and Nji are given by the relations (TT9]) and fTSOj) with n = L + 1. We 
have verified the validity of the relation ( 12"4"|) by explicit computing the partition function 
Z(X\, . . . , Xl) for L = 1,2,3,4,5 using the definition (flUjl . and in what follows we shall 
discuss the properties of the Eq. (T24|) as well as some particular solutions. 



4 The homogeneous limit 

The results of the previous sections focus on an inhomogeneous lattice whose statistical 
weights are parametrized by variables {A^} and {//&}• In this section we analyze explicitly 
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the homogeneous limit 

Afe ->■ A /i fc -)> /i (25) 

in the Eq. fl2^|) for L — 1,2 which already unveils the properties that uniquely determine 
the partition function ffTUl) . 

The case L — 1 is trivial since the partition function Z(A) can be promptly written 
down from the definitions ([2]), ([5]) and (TTOj) . However it is worthwhile to look at (124)) with 
L = 1 for illustrative purposes. In that case we have the equation 

M 1 Z(X 2 ) + M 2 Z(\i) + N 21 Z(X ) = , (26) 

and by inspecting the relations ( fT9|) and (l20i) we find in general that the limit /i^ — > fi can 
be easily obtained. By way of contrast the limit A& — > A is highly non-trivial due to the 
presence of poles in the functions Mj and Nji when two variables A& coincide. In spite of 
that, this limit can be taken using L'Hopital's rule and we shall see that the polynomial 
structure and the asymptotic behavior discussed in the appendix A also plays an important 
role in the unique determination of the the partition function Z(Xi, . . . , Xl). 

Considering the limits A ,Ai,A 2 —> X and fix — > ii in the Eq. f l2"B"j) we obtain the 
following differential equation, 



dZ x 
(q + q' 1 ) J dx 2 



2qx 



^QX 2 2 

+ q x 



(q + q 1 

expressed in terms of the variables x = e 2( - x ~^ and q = e 7 . The Eq. (127|) possess the general 
solution 

Z(x) = K l + K 2 [(l + q 2 )(xq 2 -x- 1 )-4q 2 \ogx] , (28) 

where K\ and K 2 are arbitrary constants. However, the polynomial structure discussed in 
the appendix A requires that K 2 = and we are thus left with Z(x) = K\. The constant 
K\ is then fixed by the asymptotic behavior, see appendix A, and for L = 1 we end up with 

Z{x) = (29) 

As mentioned before, the case L = 1 is a trivial case and we now turn our attention to the 
case L = 2. In that case the Eq. ffM§ reads 

M x Z(X 2 , A3) + M 2 Z(X U A 3 ) + M 3 Z(X U A 2 ) 

+ iV 21 Z(A , A 3 ) + N 31 Z(X , A 2 ) + 7V 32 Z{X , Ai) = (30) 



and again we can use L'Hopital's rule in order to evaluate the limit Ao,Ai,A2,A3 — > A and 

^1^2 f 1 - In terms of the variables Xi = e 2 ^ Xi ~^ l \ we have that Z[x\^X2) = z ^l^ where 

the function Z(x\,X2) is a polynomial of degree 1 in each variable separately as discussed in 

the appendix A. Thus, in the homogeneous limit, the function Z is a polynomial of order 2 

d n Z 

in the variable x. This fact implies in — — = for n > 2 and taking this into account, when 

dx n 

considering the homogeneous limit in (130 p . we obtain the equation 

So(z)Z(*) + + $2(*)^^ = (31) 

where 

$ ( x ) = -4g 2 (l + q 2 + q 4 ) + 6q 4 (l + q 2 )x + 12gV - 6g 6 (l + q 2 )x 3 

$i(x) = -(1 + 2q 2 + 2q A + g 6 ) + 4g 2 (l + q 2 + q 4 )x - 12q 6 x 3 + q 4 (-l + Aq 2 + Aq 4 - q 6 )x 4 
<$> 2 (x) = (1 - q 2 - q A + q 6 ) x - 2q 2 (l - 2q 2 + q 4 )(l + q 2 x 2 )x 2 + q 4 (l - q 2 - q 4 + q e )x 5 . 

(32) 

Looking to the polynomial solution of f[3"Tj) . we can generically write 

Z(x) = k 2 x 2 + kix + k (33) 



where k , k± and k 2 are arbitrary constants. In this way by replacing f l3"5|) into f[3"Tj) we find 

fc = -2 and fci = - — — — , (34) 
<T (1 + <r) 

and the constant k% is then fixed by the asymptotic behavior of Z which can be found in the 
appendix A. It turns out that 

h = ~ q~ l ?{l + q 2 ) , (35) 

lb 

and we close this section summarizing our results. 

Although we have considered only solutions in the homogeneous limit, this analysis 
suggests that the partition function Z(\\, . . . , A^) is determined by the following properties: 

(i) Functional relation 



(ii) Polynomial structure: 



Z(\ u ...,\l) = 2{X1 l '-' Xl) (36) 



n 



L-l 
2 



Xn 

8=1 

where Z{ polynomial of degree L — 1 in each variable Xi separately; 

10 



(iii) Asymptotic behavior Z(x\, . . . ,xl) ~ - ^ a ^ [L]q 2 ^ ( x i ■ ■ - x l) l 1 as Xi — > oo. 

Here [L] g 2\ corresponds to the g-factorial function whose definition is given in the appendix 
A, and in the next section we shall discuss the validity of the properties (i)-(iii) for the 
non-homogeneous case. 



5 Non-homogeneous solutions 

In order to analyze the role of the conditions (i)-(iii) given in the previous section we shall now 
consider the case with /i^ = 0, while are kept arbitrary, in addition to the homogeneous 
limit discussed previously. Let us start again looking to the case L = 1 where the partition 
function Z is constrained by the Eq. f )26|) . In that case the polynomial structure (ii) implies 
that Z is a constant and it solves the Eq. f }26]) due to the identity 



M 1 + M 2 + N 21 = , (37) 

which holds even when ^ is kept arbitrary. That constant is then determined by the property 
(iii) similarly to the homogeneous case discussed in the previous section. 

For the case L = 2, the partition function Z(Ai,A2) satisfies the Eq. (1301) and the 
polynomial structure (ii) implies in the following general form for the partition function, 

i 

Z(X 1 ,X 2 )= h m , n e mX ^ nX > (38) 

m,n=— 1 

where h m ^ n are arbitrary coefficients. By replacing the expression ( 138]) in the Eq. (1301) and 
equating the coefficients we find constraints for h m n which are solved by 

h-i-i = -y 1 h 00 = 
q A 

h-ifi = h 0; i = 

2 2 

(i + r) (i + r) 

fco,-i = />i,o = 0, (39) 

while the leading term coefficient is then determined by the condition (iii), 

h^ = ^-{ q ~ q - l f{l + q 2 ). (40) 
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The fact that the Eq. (1241) determines all coefficients of the polynomial Z(X\, . . . , Xl), 
except for one, would be expected since the Eq. ( I24p is invariant under the transformation 
Z — > Q Z where Q is an arbitrary constant. Thus the role of the condition (iii) is the 
determination of the leading term coefficient of the polynomial Z(Xi, . . . , Xl)- The partition 
function Z(Xi, A 2 ) is also expected to be symmetric in the variables Ai and A 2 , i.e. Z(X\, A 2 ) = 
Z(X 2 ,Xi), due to the definition ( ITU]) and the commutation relations (1TTI) . Though we have 
not made this assumption, the solution of the Eq. ( 130|) given by (138|) and (l39j) is manifested 
a symmetric polynomial reflecting that this property is already incorporated in the Eq. (124"]) . 

Turning our attention to the case L = 3 the Eq. (12~4"|) reads 

M 1 Z{X 2 , A 3 , A 4 ) + M 2 Z{X 1 ,X 3 ,X i ) + M 3 Z{X 1 ,X 2 ,X i ) + M 4 Z{X 1 ,X 2 ,X 3 ) 
+ N 21 Z(X , A 3 , A 4 ) + N 31 Z(X , A 2 , A 4 ) + N 41 Z(X , A 2 , A 3 ) 
+ N 32 Z(\ Q ,\ 1 ,\ 4 ) + N 42 Z(\ ,\ 1 ,\ 3 ) + N 43 Z(\ ,\ 1 ,\ 2 ) = (41) 

where the coefficients Mj and Nji are given by the Eqs. ( 1191) and (l2"Uj) . Restricting ourselves 
to the polynomial solution of ( l4ip characterized by the property (ii), the partition function 
Z(Xi, A 2 , A3) can be generically written as 

2 

Z(Ai, A 2 , A 3 ) = Yl h ^n,o e mAl+ " A2+oA3 , (42) 

m,n,o=— 2 

which can be directly substituted in the Eq. ( l4Tj) . By doing so we find constraints fixing 
all the coefficients h m ^ n ^ except for one, i.e. the leading term coefficient /i 2;2i2 , which is then 
fixed by the property (iii) 

^2,2,2 = ^(q ~ q^Hl + q 2 )(l + q 2 + q 4 ) ■ (43) 

In order to avoid an overcrowded section we have collected in the appendix B the non- 
null coefficients h m ^ nfi contained in the expression ( I42p . As we also observed for the case 
L = 2, the resulting expression (T42l is symmetric under the exchange of variables as a direct 
consequence of the Eq. (T24l) . 

We close this section remarking that although we have not rigorously proved that 
the conditions (i)-(iii) uniquely defines the partition function ([TO]) , the explicit computations 
performed in this section and in the previous one, using solely the properties (i)-(iii), strongly 
suggests that those properties completely determines the partition function Z(X±, . . . , Xl)- 
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6 Concluding remarks 



In this paper we have obtained a functional equation for the partition function of the six 
vertex model with domain wall boundary condition. The main ingredient in our derivation 
is the Yang-Baxter algebra pointing out a novel branch of exploration of this algebra. 

Although we have explicitly verified the validity of the functional equation fl24|) for 
L = 1,2,3,4,5 , it would be still interesting to provide a proof of the Izergin determinant 
representation [9]. It is also worthwhile to remark that functional equations for this partition 
function have already been presented in the literature [22] when the anisotropy parameter 
q is a root of unity. The Eq. (T2~4"l) is valid for general values of the anisotropy parameter 
q and certainly it would be interesting to investigate the connections between our equation 
and the ones of [22J. 

Different equations and representations for this partition function have also being in- 
vestigated in the literature [7J QUI [22j [23] usually obtained from Izergin determinant formula. 
Our analysis do not rely on the determinant formula and we hope the Eq. (I24p still allows 
one to investigate alternative representations. 

A careful examination of Eq. ([24"!) . or the particular case depicted in f l30l) . suggests 
that we have a relation between the partition function in the complete homogeneous limit 
and the one in the partial homogeneous limit. This partial homogeneous limit has been 
discussed in [23] associated to the proof of the refined alternating sign matrix conjecture, 
and in [25] related to the six vertex model artic curve. The examination of (124"|) then suggests 
that such an equation relating partial and complete homogeneous limit could be obtained 
by considering Ai, A 2 , . . . , Xl+i — > A while keeping A fixed. 

Within the framework of the Quantum Inverse Scattering Method, the study of par- 
tition functions with domain wall boundary conditions shares many aspects with the com- 
putation of scalar products and correlation functions [13]. Given this similarity we hope 
to report on functional equations for scalar products and correlation functions in a future 
publication. 
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Appendix A: Polynomial structure and asymptotic be- 
havior. 



As it was shown in [I] and [15], the definitions (J2]), (EJ) and (jHJ) allow us to study the 
dependence of the operator B(X, {/ifc}) with the spectral parameters. In terms of the variables 
%i = e 2 ( Al ~^, it turns out that 

1 



B(Xi,{(Mk}) 



x. 



0, 



r(0 



(a.i; 



with L operator coefficients f a . The leading term fi-i can be written down explicitly due 
to the structure of (J2J) and and it is given by 



= 2- V^V - 1) exp 

where the operators Pj are defined as 

j'-i 



i=i 



K ® X- <8) (g) K- 1 . 



Together with X~ 



1 




the matrices 



(A.2) 



(A.3) 



X~ 





1 o 



and K 



i 



generate the g-deformed su{2) algebra, 

KX ± K~ l 

[x\x-] 



q ± X ± 
K 2 - K- 2 
q-q- 1 



(A.4) 



(A.5) 



Looking to the definition ([TTJ]) we can readly see that the dependence of Z(Xi, . . . , Xl) 

l 

with the spectral parameters follows directly from the dependence of Y\ B(Xi, {fik}) with 



i=l 



{Xj} and {fj.k} since the vectors |0) and |0) do not depend on them. Thus considering the 
definition (flUj) and (1A.1I) we can write 



Z{Ai, . . . , A L ) = 



(A.6) 



n 

i=i 
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where Z(xi, . . . , xl) consists of a polynomial of degree L—l in each variable Xi independently. 

In their turn the coefficients f^-i govern the asymptotic behavior of Z. In the limit 
X{ — > oo only the leading term coefficients contribute and from (lA.lj) and (1A.6I) we obtain 

L 

Z(x u ...,x L )^(0\Y[f^ 1 \0) {x l ...x L ) L ~ 1 as x^oo. (A.7) 

i=l 

L 

The product J^J can be worked out using the following properties exhibited by the 



operators Pj, 



i=l 



PiPj = q 2 P ] P l (i<j) (A.f 

>2 



Pf = 0, (A.9) 

which can be derived with the help of the relations (lA.3j) - f)A.5j) . 
Considering the property (1A.9[) . we then have 

f[ fS-i = 2- L2 q ^ (q 2 ~ If Y, P«P« -Pa, (A.10) 

i=l { aj } 

where denotes a summation over {cij} with each index aj ranging from 1 to L under 

{ aj } 

the constraint a\ ^ a 2 ^ ■ ■ ■ ^ a^. The terms of P ai P a2 ■ ■ ■ Pa L can be ordered using the 

{%} 

relation (IA.8[) in order to achieve a common element. Thus we find 

Y P ai P a2 ...P aL = (l + q~ 2 ){l + q- 2 + q- 4 ) ■ ■ ■ (1 + q' 2 + q~ 4 + ■ ■ ■ + q~^ X) )PxP% . . . P L . 

{%} 

(A.11) 

Now in order to derive an explicit expression for the leading term coefficient in (|A.7[) , we only 
need to compute (0| P\Pi . . . Pl |0). This task can be directly performed since P\P 2 . . . Pl 
consists of a tensor product of local operators. Using (I12p and (1A.4j) we thus obtain 

(0| P 1 P 2 ... P L |0) = . (A.12) 



Gathering our results so far, in particular the relations OA.lOp . flA.lip and (IA. 12[) . we 
are left with 

(0\f[ftx\0)= iq ~ 2 % 1)L [L%> (A.13) 

i=i 

where [L\] q2 denotes the q- factorial function defined as 

[L] q2 \ = 1(1 + q 2 )(l + q 2 + g 4 ) . . . (1 + q 2 + ■ ■ ■ + q 2 ^) . (A.14) 
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Appendix B: The coefficients h m ^ 

In the table below we present the non-null coefficients fa m , ni0 of the polynomial (Pf2]) solving 
the Eq. flUJ). 



(m, n, 6) 


f^m ,n,o 


(m, n, o) 


^2,2,2 


(-2,-2,-2) 




(0,0,0) 


Cl — 8(7 2 — 34*o 4 — 8(7 6 -l"r7 8- l 


g 4 (l+ g 2 )(l-q+g 2 )(l+ g + (? 2 ) 


(-2,-2,0) 


— 3(l+g 2 ) 


(0,0,2) 




9 4 (l-q+9 2 )(l+g+g 2 ) 


(l-9+g 2 )(l+g+? 2 ) 


(-2,-2,2) 


q 
O 


(0,2,-2) 


-(l+10g 2 +g 4 ) 


g 2 (l-g+ g 2 )(l+g+g 2 ) 


g 2 (l+g 2 )(l- g + g 2 )(l+g+<j 2 ) 


(-2,0,-2) 


-3(l+g 2 ) 


(0,2,0) 




9 4 (l-g+9 2 )(l+g+g 2 ) 


(l- g +g.2)(l +g+ y2) 


(-2,0,0) 


1 


(0,2,2) 


-3(l+g 2 ) 


g 2 (l-g+g 2 )(l+ g +q2) 


(l- g + g 2)(l +g+ y2) 


(-2,0,2) 


-(l+10a 2 +a 4 1 


(2,-2,-2) 


Q 
O 


g 2 (l+9 2 )(l-<?+g 2 )(l+g+<? 2 ) 


<? 2 (l-g+g 2 )(l+<?+g 2 ) 


(-2,2,-2) 


3 


(2,-2,0) 


-(l+10 g 2 +g 4 ) 


g 2 (l-g+ g 2 )(l+g+q 2 ) 


g 2 (l+g 2 )(l- g +g 2 )(l+g+<j 2 ) 


(-2,2,0) 


-(l+10g 2 +<? 4 ) 


(2,-2,2) 


3 




(l- 9+g 2)( 1+9+? 2) 


(-2,2,2) 


3 


(2,0,-2) 


-(l+10g 2 +g 4 ) 


(l-<Z+<Z 2 )(l+?+<7 2 ) 


g 2 (1+g 2 )(1 _ g+(/ 2 )(1+g+9 2 ) 


(0, -2, -2) 


-3(l+g 2 ) 


(2,0,0) 


12 




(l- 9+ g2)( 1+9+? 2) 


(0,-2,0) 


12 


(2,0,2) 


-3(l+9 2 ) 


g 2 (l-g+ g 2 )(l+g+<? 2 ) 


(l- 9+ g2)( 1+9+? 2) 


(0,-2,2) 


-(l+10g 2 + 9 4 ) 


(2,2,-2) 


3 


q 2 (l+ g 2 )(l-q+g 2 )(l+g+g 2 ) 


(l- g +g 2 )(l+g+g 2 ) 


(0,0,-2) 


12 


(2,2,0) 


-3(l+g 2 ) 


9 2 (l-g+? 2 )(l+9+3 2 ) 


(l- g +g 2 )(l+ g +g 2 ) 
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